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We implement the effects of disorder on a holographic superconductor by introducing a random chemical 
potential on the boundary. We demonstrate explicitly that increasing disorder leads to the formation of is¬ 
lands where the superconducting order is enhanced and subsequently to the transition to a metal. We study the 
behavior of the superfluid density and of the conductivity as a function of the strength of disorder. We find 
explanations for various marked features in the conductivities in terms of hydrodynamic quasinormal modes 
of the holographic superconductors. These identifications plus a particular disorder-dependent spectral weight 
shift in the conductivity point to a signature of the Higgs mode in the context of disordered holographic super¬ 
conductors. We observe that the behavior of the order parameter close to the transition is not mean-field type as 
in the clean case, rather we find robust agreement with exp(—A \T — with v = 1.03 ± 0.02 for this 

disorder-driven smeared transition. 


Introduction: The suppression of conductivity due to dis¬ 
order, known as Anderson localization, is one of the most 
striking transport phenomena that involves quantum behav¬ 
ior m. The application of localization ideas to superconduc¬ 
tors was for a long time governed by Anderson’s theorem stat¬ 
ing that nonmagnetic impurities have no significant effect on 
the superconducting transition since Cooper pairs are formed 
from time reversed eigenstates, which included disorder 12 . 
Anderson’s idea applies only to weakly disordered systems 
(with extended electronic states) and weak interactions; the 
role of interactions has been further discussed by Ma and Lee 
in IJI. The mechanism put forward in 13 states that strong 
disorder gives rise to spatial fluctuations of the order parame¬ 
ter along with its suppression in comparison to its value in the 
clean system. The existence of spatial fluctuations of the order 
parameter has recently been corroborated and argued to be the 
central mechanism in the superconductor-insulator transition 
liiia. Various experiments have indicated that high-Tc su¬ 
perconductors are intrinsically disordered 012 emphasizing 
the need for a better understanding of the interplay between 
superconductivity, strong interactions and disorder. In this 
manuscript we tackle this interplay using holographic meth¬ 
ods. 

The Anti de Sitter/Conformal Field Theory (AdS/CFT) 
correspondence has already succeeded in constructing holo¬ 
graphic versions of superconductors 0 |3 (for reviews see 
inaini). Some holographic models have been argued to be 
relevant to describe materials like the cuprates ifT^ [131 since 
their dynamics is believed to be largely governed by their 
proximity to a quantum critical point 1141 . 

In previous work 1131161 we investigated mild disorder 
in holographic superconductors and found an enhancement 
of the superconducting order parameter. Other recent works 
studying disorder holographically include inHisi. Here we 
explore the regime of large disorder and provide a detailed de¬ 
scription of the superconductor-metal transition. Although we 
work in a setup where the normal phase describes a metal, it is 


also possible to consider insulating setups through AdS/CFT 

EoMl. 

In this paper we report various properties of the super¬ 
conducting transition including the averaged and spatially re¬ 
solved behavior of the condensate, superfluid density and con¬ 
ductivity and provide an explanation in terms of the rele¬ 
vant hydrodynamic modes. We observe that the disorder- 
driven transition is smeared and the averaged order param¬ 
eter (here the average (•) is over the spatial coordinate or, 
equivalently, disorder realizations) vanishes very quickly as 
(O) ~ exp {-A\T - we And i/ = 1.03 ± 0.02. We 

expect that our explicit result for this exponent will stimulate 
alternative approaches to disorder-driven transitions to pro¬ 
vide a quantitative characterization of the transition. 

Disordered holographic superconductor: To build a 
noisy holographic s-wave superconductor in 2 h- 1 dimensions 
we consider, following 0 , the dynamics of a Maxwell held 
and a charged scalar in a fixed metric background: 

The system is studied on the Schwarzschild-Ac/S '4 metric: 

^ + dx'^ + dy^^ , 

f{z) = (1) 

where we have set the radius of AdS, R = 1, and the horizon 
at z/j = 1. This system is dual to a 2 h- 1 CFT living on the 
boundary of AdS^, and the C/(l) gauge held realizes a con¬ 
served current. The temperature of the black hole is identified 
with that of the held theory, and by fixing the horizon radius 
we are making use of the rescaling symmetry of our theory to 
work in units of temperature. We take the following (consis¬ 
tent) Ansatz for the matter fields: 

T'(a:, z) = tpix, z), A = (f){x, z) dt, 


( 2 ) 
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where z) G K. The resulting equations of motion read 

dU+jdU-^</> = 0 , ( 3 ) 

dli^ + 7 + ("7 “ 7) + 72 - ”^) V’ = 0. 

(4) 

In what follows we choose the scalar mass m? = —2, corre¬ 
sponding to a dual operator of conformal dimension A = 2. 

The main rule for reading the AdS/CFT dictionary states 
that field theory information is extracted from the boundary 
values of the gravity fields. The UV (z — 0) asymptotics of 
Eqs. ffl lead to 

4>{x, z) = fj.{x) + p{x) z + (x) z^ + o(z^), (5) 

^(x, z) = 4^^\x) z + (x) ^ + o{z^), (6) 

where [i{x) and p{x) correspond, in the dual field theory, 
to space-dependent chemical potential and charge density re¬ 
spectively. The functions 4^^\x) and tp^^Hx) xie identified, 
under the duality, with the source and VEV of an operator of 
dimension 2. Imposing pj^^Kx) = 0 in Eq. (j^ corresponds 
to spontaneous breaking of the C/(l) symmetry with order pa¬ 
rameter O oc In the IR {z ^ Zh = 1) regularity 

implies that At vanishes at the horizon. 

To mimic the choice of random on-site potential used orig¬ 
inally by Anderson in HI we implement disorder by introduc¬ 
ing a noisy chemical potential; 

p{x) = f^o + ^u> X! cos{kx + Sk), ( 7 ) 

k—ko 

where 6k is a random phase for each k, and w is a free parame¬ 
ter that determines the strength of the disorder. Notice that the 
system is homogeneous along the remaining spacelike direc¬ 
tion y. We discretize the space, and impose periodic boundary 
conditions in the x direction, leading to k with values; 

Oyr Ic 

kn = — (n + 1) with 0 < n < N = , (8) 

L ko 

where L is the length in the x direction of our cylindrical 
space. Our noise is a truncated version of Gaussian white 
noise where the highest wave number A:* takes the role of the 
inverse of the correlation length for the chemical potential. 
More details on the properties of this choice of disorder can 
be found in ifT^ . 

Summarizing, we think of fcg as the inverse system size and 
of kf, as the inverse correlation length, and we work in the 
regime k^/T 1. More precisely, for most of the simula¬ 
tions we take L = 207r, and fc* = 1, but we checked the 
stability of the results for lengths up to L = SOtt. 

To find solutions describing disordered superconducting 
states we integrate the equations of motion (|^[^ for different 
realizations of disorder characterized by sets of the random 
phases 5k- We define the expectation values for the different 


observables by averaging over different realizations. Typical 
plots correspond to about fifty realizations. 

The order parameter and superconducting islands: Let 
us focus on the behavior of the order parameter as a function 
of the chemical potential p and the dimensionless strength of 
disorder w. We will pay special attention to the minimum of 
the condensate Omin throughout the sample which we expect 
to be related with the DC conductivity. 

InlBKIl it was shown that noise enhances the spatial av¬ 
erage of the order parameter. However, already Ma and Lee 
found in 13 that strong disorder gives rise to spatial fluctu¬ 
ations of the order parameter along with its suppression in 
comparison with its value in the clean/homogeneous system. 
In Eig.[T|we plot the spatial dependence of the order parameter 
at a temperature T — 0.8177^°, where 77 ^° is the homo¬ 
geneous critical temperature. We find that for strong enough 
noise regions akin to islands appear in the system. 



Figure 1; The condensate as a function of the coordinate x for 
three realizations with strength w = 0.1,0.9,3, and the same set 
6k (orange, yellow and purple lines respectively), at a temperature 
T = 0.81T“^°. This plot shows the appearance of islands, that is, 
of spatial fluctuations in the condensate. 

The enhancement of superconducting properties reported 
in ifTSl ri6l has some precedent in the condensed matter lit¬ 
erature. The role of superconducting islands was already ex¬ 
plicitly mentioned in 12^ . and further analyzed in Q. In that 
work it was shown that for each realization of disorder, there 
are spatial regions where the local upper critical field exceeds 
the system-wide average value. These regions form super¬ 
conducting islands weakly coupled via the Josephson effect. 
At low temperatures, proximity coupling is long ranged and, 
thus, global superconductivity may be established in the sys¬ 
tem. Similar arguments were also advanced in 0. We see 
that precisely this mechanism seems to be at play in our holo¬ 
graphic model. 

To characterize the possible appearance of islands we study 
the minimum of the condensate O^nin for different values of 
the parameters that characterize the system. In Eig.j^we plot 
Omin (dashed line) as function of the disorder strength w at 
constant temperature T = 0.8177^°. This figure shows that 
Omin is a holographic version of the original suggestion of 
Ma and Lee 0 that strong disorder gives rise to spatial fluc¬ 
tuations of the order parameter along with its suppression in 
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Figure 2: Superfluid density, Us (solid line) and minimum of the 
condensate, Omin (dashed line), as a function of disorder, w, for 
= 0.81. The value of Omin has been normalized to unity 

at ui = 0. 

comparison with its value in the clean/homogeneous system. 
Moreover, for large values of the disorder strength we observe 
an exponential tail for Omin- This means that we have some 
leaking between islands of superfluid leading to a finite value 
of the condensate in between. We relate this to the fact that 
AdS suppresses large momentum contributions to the order 
parameter as reported in ifTSl [Thll . 

Conductivities: The holographic computation of conduc¬ 
tivities requires the study of fluctuations on top of the back¬ 
ground. Here we focus on the conductivity of the C/(l) bro¬ 
ken superconducting phase, where the perturbations couple to 
the noise due to the nontrivial scalar held. Given the sym¬ 
metries of the setup, there are two different conductivities one 
can study; the one in the direction parallel to the noise and that 
orthogonal to it. The orthogonal conductivity is insensitive to 
the disorder, and we thus focus on the parallel one. 

To compute the electric conductivity we consider homo¬ 
geneous perturbations of the gauge held oscillating with fre¬ 
quency oj: 

(9) 

and on the boundary (z = 0) we impose fti{x) = 1 (where 
is the held strength of a^, and i runs over the spatial direc¬ 
tions X, y), thus sourcing a constant electric held in the dual 
held theory. Then, following AdS/CFT, the conductivity is 
computed as 


held x(a:, z) + iri{x, z), and we must then solve the four cou¬ 
pled linear equations of motion ( [T5] l. Interestingly, as we will 
show in Fig. the fields (%, p) show up as hydrodynamic 
modes in the computation of conductivities. 

When looking for the low frequency behavior of the con¬ 
ductivity we And that, as in the homogeneous phase, we have 
a pole in the imaginary part of cr. This translates, through 
Kramers-Kronig rules, into a (numerically invisible) delta 
function in the real part, giving 

a Ki Us + ... , ( 11 ) 

where Ug is the superfluid density. Notice that for axx, a near¬ 
boundary analysis of the equations of motion shows that, due 
to current conservation, the DC conductivity must be homo¬ 
geneous, and therefore rig is a constant independent of x. 

The dependence of the superfluid density on the strength 
of disorder is shown in Fig. together with the evolution of 
the minimum value of the condensate. For low w, we verify 
that Us does not change much in total agreement with Ander¬ 
son’s theorem even if we are describing a strongly coupled 
system. The persistence of the transport properties for weak 
disorder, even in strongly coupled systems, has been argued 
in the metal-insulator transition in Il24l . Fig.j^shows that for 
large w, both the superfluid density and the minimum of the 
condensate decay exponentially. This illustrates the disorder- 
driven transition to the normal (metallic) phase. 



Figure 3: Real part of the AC conductivity for r/T™^° = 0.45. The 
black line corresponds to the homogeneous case, whereas the purple 
and yellow lines denote the disordered w = 1 and tu = 2.4 cases 
respectively. The dimensionless frequency Cj is proportional to lj/T. 


ctu{x) 




idzai{x, z) 
uja^{x,z) ^^0 


( 10 ) 


where there is no summation over i, and we impose ingoing 
boundary conditions for ai(x, z) at the horizon. 

For the conductivity in the direction of the noise we must 
consider the linearized equations for the spatial component 
of the gauge held ax{x, z). Unfortunately, the x dependence 
couples Ux to the perturbation of the temporal component of 
the gauge held ao(x, z) and to the perturbations of the scalar 


We present the real part of the AC conductivity in Fig. 
comparing the homogeneous case with the results for two dif¬ 
ferent values of the strength of disorder. The two key fea¬ 
tures absent in the homogeneous case are; the presence of 
resonances for small frequencies and a shift of the spectral 
weight at large frequencies. The small frequency resonances 
are a direct result of coupling of new fluctuating fields for 
an cc-dependent chemical potential. These resonances can be 
understood as related to the holographic quasinormal modes 
studied in ll^ l26l and shown to affect the conductivity at 
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nonzero superfluid velocity in lIZTll . The shift in the spec¬ 
tral weight clearly depends on the strength of disorder, and 
we view it as evidence of the Higgs mode associated with the 
spontaneous breaking of the C/(l) symmetry ll^ . The Higgs 
mode corresponds to oscillations of the modulus of the order 
parameter, and it is predicted ll29l to cause an excess of the 
electrical conductivity at sub-gap frequencies. In clean BCS 
superconductors the mass of the Higgs mode results in a gap 
of its contribution to the conductivity which is of the order of 
the BCS gap, and thus makes it difficult to detect the Higgs 
mode through its effect on the conductivity. However, it was 
shown in ll^ that disorder suppresses that gap, giving rise to 
an observable excess of the AC conductivity at sub-gap fre¬ 
quencies with respect to the standard BCS prediction (based 
on the measurement of the energy gap of Cooper pairs) 1281 . 
The excess observed in l28l is qualitatively similar to what we 
observe in Fig.|^at frequencies 2 < w < 4 where the conduc¬ 
tivity of the disordered system (yellow line) is higher than that 
of the clean one (black line). 

To verify that the small frequency resonances in the con¬ 
ductivity correspond to the gapless modes studied in l2^ . we 
compute the effective velocity, Vg, of the first of these reso¬ 
nances assuming the dispersion relation 

oj = Vs{T,w) k , (12) 

and taking k to be equal to fco, where fcg = 27r/L is the small¬ 
est wave number in the sum Q. We then compare the evolu¬ 
tion of Vs with the temperature with the findings of l26l . The 
result is plotted in Fig. which shows that Vb(T,w) follows 
the standard temperature dependence discussed in 1^ |26l 
(corresponding to the solid black line in the graphic). What 
is new in our case is its dependence on the disorder. For small 
temperatures the role of disorder is suppressed. For higher 
temperatures we see that Vg decreases with increasing disor¬ 
der strength except very near the critical temperature, Fig.|^ 



’H '=0 


Figure 4: Effective velocity Vs {T, w) corresponding to the first res¬ 
onance of the conductivity for w = 0.1 (blue dots) and ui = 1 (red 
dots). The black solid line presents the result of 1261 for the speed of 
sound of the sound mode in holographic superconductors. 

Smeared phase transition: Finally, let us turn to the anal¬ 
ysis of one of the most universal properties in phase transi¬ 
tions: the behavior of the order parameter close to the transi¬ 


tion M- As a first approach to disordered phenomena, the 
fate of a particular clean critical point under the influence of 
impurities is controlled by the Harris criterion lf32l that gen¬ 
eralizes the standard power-counting criterion to random cou¬ 
plings. As explained in ll33l our chemical potential disordered 
along one dimension, and thus perfectly correlated along the 
remaining spatial direction, introduces relevant disorder. This 
places our setup in the class of systems where quenched dis¬ 
order can lead to exotic critical points where the conventional 
power-law scaling does not hold llMl - l3^ . Moreover, in some 
cases disorder has been shown to cause the formation of rare 
regions undergoing a phase transition independently from the 
rest of the system. In terms of the average of the order param¬ 
eter this results in a smeared phase transition characterized 
by an exponential scaling l|37H39l. One would expect the is¬ 
lands of conductivity that form in our system to play the role 
of these rare regions, and consequently to smear out the phase 
transition. 



Figure 5: The disorder-driven transition is smeared. Plot of the spa¬ 
tial average of the condensate close to the critical temperature for 
w = 0, 0.1, 0,2 and 0.4 (black, purple, orange, and yellow lines, 
respectively). 

Fig. shows the transition for the clean, homogeneous 
case 1917 and for three disordered cases (w — 0.1, 0.2,0.4). 
We confirm that for the disordered transition the order pa¬ 
rameter behaves as (O) ^ exp {—A\T — Tc\~^) in Fig. 
Notice that Tc is different for each value of w, and is al¬ 
ways higher than By taking the logarithmic derivative 

{O'[T))/{0{T)), this figure shows linear fits that allow us 
to determine the exponent v, obtaining v = 1.03 ± 0.02 in¬ 
dependent of the disorder w. Explicitly, the fits in Fig. are 
y = —2.26 — 2.02a; for w = 0.4, y = —3.13 — 2.01a; for 
w = 0.2, and y = —4.21 — 2.06a; for w = 0.1, and notice 
that the slope of these fits corresponds to —{v + 1). The fits 
also show that the coefficient A of the expression above in¬ 
creases with the disorder strength w, in agreement with the 
optimal fluctuation theory predictions ll37l[39l . In particular, 
after assuming « 1, these fits result in the following values 
of A\ 0.01 at ic = 0.1, 0.04 at m = 0.2, and 0.1 at w = 0.4. 
The ranges of the fits, shown by solid black lines in Fig.[^ are 
determined on one side by the region where the condensate 
behaves as in the homogeneous case which as Fig. shows is 
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for temperatures around .95Tc. On the other side the bound 
arises from the numerical instability introduced by exponen¬ 
tially small values of the condensate. Notice that the points at 
temperatures above are rare, and many realizations are 

needed in order to have enough statistics to capture their be¬ 
havior. Moreover, those points correspond to very low values 
of the condensate, and small errors in {0{T)) affect strongly 
the value of {O'[T))/{0{T)). Finally, one should notice that 
the UV cutoff A:* in our disordered chemical potential bounds 
the disorder distribution. Consequently, one expects to hnd, 
as we do, a finite disordered critical temperature T“ above 
which only the normal phase exists Ell 139). 



Log(r,-T) 

Figure 6: The disorder-driven transition is smeared. We plot 
log versus log(rc — T) for w = 0.1, 0.2, 0.4 (purple, or¬ 

ange, and yellow respectively). The solid black lines denote linear 
fits to the data, and the vertical dashed lines indicate the critical tem¬ 
perature of the homogeneous case 

Conclusions 

Superconducting islands in holography: We argued that the 
formation of islands in the context of holographic disorder- 
driven superconductors is, as in experimental and numerical 
analyses, the central mechanism at play in the transition. This 
behavior is responsible for the enhancement of superconduc¬ 
tivity reported in MM- The fact that the islands are coupled 
via the Josephson effect, as suggested originally in IH, causes 
the exponential decay in the tail of the superfluid density as 
shown in Fig. 

Disorder in strongly interacting systems: The interplay be¬ 
tween disorder and interactions has been the subject of vigor¬ 
ous studies in the condensed matter literature initiated most re¬ 
cently by the work of Basko et al. Il24l . In the series of works 
generated by ll24l . (see for example 11312) it is shown 
that for sufficiently strong disorder, weak interactions do not 
change the nature of the localized phase Il43ll44l . The present 
work is in the complementary region, where interactions are 
strong to begin with, and disorder is increased. We have es¬ 
tablished that weak disorder does not destroy the holographic 
superconducting state, in particular. Fig. shows that the su¬ 
perfluid density is largely unaffected for small disorder. In¬ 
stead, for strong disorder the formation of islands suppresses 
the superfluid density driving the system to the normal phase. 

Optical conductivity: We have studied the AC conductivity 


for disordered holographic superconductors, identifying the 
new low-frequency resonances with quasinormal modes of the 
holographic superconductor. Moreover, the conductivity dis¬ 
plays a disorder-dependent shift of the spectral weight highly 
suggestive of a massive Higgs excitation. Strong experimental 
evidence in favor of a Higgs mode in strongly disordered su¬ 
perconductors close to the quantum phase transition has been 
recently reported in 12 ^ . 

The disorder-driven transition is smeared: We have shown 
that under the influence of disorder the superconductor-metal 
transition changes from power law, mean field to a smeared 
transition of the form exp(—A |T — Tc\~''), with ly = 1.03 ± 
0.02. We hope that our results motivate other approaches to a 
quantitative description of this transition. 

Future directions : Understanding the properties of con¬ 
ductivities in the language of modes helps us compare with 
more traditional condensed matter methods directly. Building 
a holographic disordered superconducting thin film will help 
us tackle the Higgs mode in the strongly disordered supercon¬ 
ductor close to the quantum phase transition 1^ . Ultimately, 
due to the universality of critical exponents EH, we hope that 
our holographic approach provides a direct path to disorder- 
driven critical exponents. 
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Supporting material 

Numerical stability: Given the aleatory nature of our cal¬ 
culation, in this appendix we show the stability of the result 
of the value of i/ against the number of realizations. Namely, 
Fig. [7] shows that the result stabilizes as we increase the num¬ 
ber of realizations of disorder, i.e. the number of times we 
generate a series of random phases, 5k in Eq. 0. Most of our 
simulations for determining v involve over 100 realizations. 

Homogeneous versus smeared condensate: In Fig. we 
show that the homogeneous transition follows a power law. 
We present it in the same plot as a smeared transition to il¬ 
lustrate how our numerical analysis distinguishes between the 
homogeneous and the disordered transition. Note that for an 
homogeneous condensate Oh = Ah\T — Tc\^, with /? = 1/2 
El, hence one can write 

log (^) =-log |T-Tc|-F cons, (13) 

while for the disordered case O = Ag exp(—A|T — Tc\~'^), 
and one has 

log = -{v + 1) log |T - Tel -F cons . (14) 
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Equations of the perturbations: The equations of motion 
for the perturbations relevant for computing the electric con¬ 
ductivity read: 


-4.0 -3.5 -3.0 -2.5 -2.0 

Log(r,-T) 

Figure 7: The condensates for w = 0.4 with four different numbers 
of realizations: purple corresponds to 10 realizations, orange to 25, 
yellow to 75, and black to 150. Notice that except for the leftmost, 
high temperature, data points, the value of the condensate stabilizes 
for a number of realizations above 25. 




Figure 8: The homogeneous condensate (black dots) and the smeared 
condensate for w = 0.1 (purple). The red, and blue lines correspond 
to the linear fits y = —0.71 — 1.00®, and y — —4.21 — 2.06® 
respectively. 


2{iujr] + aoi} + 2(l>x)ip - z'^ {f d^^ao + ioJ = 0 , (15a) 

2 

y flj, -F f f dzax + dzzaz; - iujdz^ag) - 2 - ip dxV -F d^ipn) = 0 , (15b) 

-2{ioj'n + ipao)(p-z'^(p'^X+ {-z^ f) x - z f (-2 f dzX + z f dzX + zf dzzX + z d^^z^x) = 0 , (15c) 


[ 2 ^ -F pP') — mP /] Tj — 2x dzy + zP ^—iujag ip — 2i ujcp 

subject to the constraint: 

-iz^ bj dzGo + f [2 Ip dzT] -2rj d^ip - d^zO-x) = 0 . 

( 16 ) 
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+ / if' dzT) -F / dzzT] - Ip dxax -2ax d^ip + d^xv)] = 0 , (15d) 
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